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Abstract 

We explore superconducting instabilities of black branes in 50(6) gauged super- 
gravity at finite temperature and finite R-charge densities. We compute the critical 
temperatures for homogeneous neutral and superconducting instabilities in a trunca- 
tion of 20 scalars and 15 gauge fields as a function of the chemical potentials conjugate 
to the three U(l) charges in 50(6). We find that despite the imbalance provided by 
multiple chemical potentials there is always at least one superconducting black brane 
branch, emerging at a temperature where the normal phase is locally thermodynami- 
cally stable. We emphasise that the three-equal charge solution, Reissner-Nordstrom, 
is subdominant to a thermodynamically unstable black brane at sufficiently low tem- 
peratures — a feature which is hidden in an equal charge truncation. 



1 Introduction 



In the holographic approach to condensed matter physics there is great freedom in how a bulk 
theory is constructed. This has the advantage of allowing the construction of holographic 
models which capture a wide variety of physical phenomena. The drawback, of course, is 
that the detailed phenomenology of such constructions typically depends on the particular 
model employed. For this reason a desirable feature of any phenomenological holographic 
programme is the existence of some 'universal result', at least for a particular class of bulk 
theories. A celebrated example is the ratio of shear viscosity to entropy density for strongly 
coupled, isotropic holographic plasmas arising from two derivative bulk Lagrangians [1, 2]. 

Holographic superconductivity on the other hand has no simple universal sector; in the 
minimal example of a single £7(l)-charged bulk scalar field there is freedom to specify the 
potential and gauge field coupling. In a bottom-up approach, whilst the small-field expan- 
sion of these functions can be easily parameterised, the behaviour away from the critical 
temperature depends strongly on their nonlinear characteristics. Turning to a top-down 
approach one may hope to overcome such difficulties. However, as observed in [3], different 
consistent sub-truncations of a particular top-down model exhibit radically different results 
for the behaviour of the superconducting order parameter with temperature. To understand 
if there are any generic features for the landscape of consistent holographic superconductors, 
one must delve deeper into larger families of consistent truncations. See [4, 3, 5, 6] for some 
interesting work along these lines. 

In this paper we take a step in this direction by considering a truncation of M = 8 50(6) 
gauged supergravity in 5D containing 20 scalars and 15 gauge fields. The truncation contains 
a three-charge family of black branes [7], dual to a normal phase, about which we look for 
the emergence of new branches of superconductors. Such branches are indicated by the 
existence of marginal modes for the various charged scalar fields at finite temperature. This 
calculation can be viewed as a first step towards understanding the sequence of preferred 
phases of the full supergravity theory. 

When cooled sufficiently the black branes we employ to model the normal phase become 
thermodynamically unstable. 1 Note that in many cases the stable black brane simply ceases 

: At high temperatures in the regime of local thermodynamic stability there are indications that this 
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to exist at low temperatures, and so the system is largely immune to any T = test 
of stability, such as Breitenlohner-Freedman (BF) bound violation for any near horizon 
AdS2 x M. 3 . We focus on superconducting instabilities that occur in the normal phase at 
sufficiently high temperatures where it is locally thermodynamically stable. We find that a 
superconducting black brane branch always exists in this region of the normal phase. 

In addition, the black brane solutions we employ for the normal phase have three R- 
charges, which allows us to explore the dependence of the various superconducting phase 
transitions on the conjugate chemical potentials. Some modes in the truncation are charged 
under a single £7(1) and do not couple to the other gauge fields at the linearised level. 
As such, turning on an additional chemical potential for another U(l) provides an effective 
holographic description of a state with a density imbalance. An example in condensed matter 
physics is an imbalanced Fermi mixture in which there is a mismatch of chemical potentials 
for spin-up and spin-down constituents [11]. Phenomenological holographic models of such 
mixtures were investigated in bottom-up bulk and probe brane constructions in [12, 13], 
where the additional U(l) plays the role of up-down spin imbalance. One feature of these 
models is that the superconducting phase ceases to exist if the imbalance is too great. 

Here we have an imbalance as a natural consequence of the three chemical potentials. 
As the imbalance for a given superconducting mode is increased, its critical temperature 
decreases until it reaches the threshold of thermodynamic instability, at finite temperature. 
However, increasing the imbalance for one mode decreases the imbalance for another, and 
is a consequence of the SO (6) symmetry. Surprisingly there is no set of chemical potentials 
for which the imbalance is simultaneously too great for all modes. 

The paper proceeds as follows. In section 2 we present the consistent truncation and 
the three-chemical potential black brane solution we will be employing as the background 
normal phase. In section 3 we calculate the field equations for homogeneous fluctuations 
of the other fields in the truncation about this background that do not source gravitational 

system is only metastable [8, 9, 10]. With the interpretation of the three R-charges as angular momenta of 
D3-branes spinning in the S 5 of the full 10D theory, this instability corresponds to the one-by-one fragmen- 
tation of D3-branes from the rotating stack. The superconducting branches we seek in this work are local 
instabilities of the black brane solutions, and whilst they may be subject to arguments vis-a-vis metastability, 
a study of the critical temperatures and chemical potential dependence is not directly affected by it. The 
family of black branes used here represents a convenient local minimum for our purposes. 
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fluctuations. The structure of the equations for the fluctuations and the instabilities found 
naturally leads to two different types of mode: modes that are continuously connected to 
charged modes are discussed in 4, while those that are continuously connected to neutral 
modes are discussed in section 5. In both cases we map out the regions in the space of 
chemical potentials which are unstable to each mode and detail their critical temperatures. 
We conclude in section 6. 

2 Truncation and the normal phase 

In addition to gravity, the bosonic sector of M = 8 50(6) gauged supergravity in 5D [14, 15, 
16] contains 42 scalars, 15 5*0(6) gauge fields as well as twelve 2-form gauge potentials. The 
scalars are organised into irreps of 50(6) which have masses m 2 L 2 = —4, —3, in an AdS$ 
vacuum, the lightest of which saturates the BF-bound. We will employ a further consistent 
truncation of this theory down to 15 gauge fields and 20 of the BF-bound saturating scalars. 
In 10D ingredients this truncation corresponds to retaining only the metric and 5-form [17]. 

This truncation provides a convenient starting point in our search for superconducting 
black brane branches. First, the truncation contains a three-parameter family of black brane 
backgrounds which we may employ as the normal phase. Second, only the lightest scalars 
remain and these should display a greater susceptibility to forming condensates in the black 
brane backgrounds. 2 Finally, we are able to make contact with special cases studied elsewhere 
at the nonlinear level. 

The 20 scalars and 15 gauge fields can be packaged as a symmetric unimodular 50(6) 
tensor Tij and an antisymmetric 50(6) tensor A^, respectively. The equations of motion for 
this consistent truncation may be derived from the following action, with the adjoint 50(6) 
indices written in matrix form: 

5 = — ^- j d b x^-g (r - ^tiT' 1 (D M T) T' 1 (£>T) + l -ti T~ l F^T' 1 - V^j (1) 

We have omitted a Chern-Simons term that will not be important for the electric, homoge- 
2 Superconducting branches for special cases of the other scalar field sectors are investigated in [5]. 
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neous superconductors sought here. The scalar potential is given by 

V = & (2trT 2 - (trT) 2 ) , (2) 

with the 50(6) gauge-covariant derivative D^T = V^T + g c [A^,T] and field strength = 
d^Ay — dyA^ + gc [A^, A„]. We set the dimensionful gauge coupling g c = 1 throughout, which 
fixes the AdS length. 

When written in 2x2 blocks, the background three-charge asymptotically AdS black 
brane solution [7] takes the form 

T = T = diag(X 1 l,X 2 l,X 3 l) (3) 
A = A dt ee diag(v4 1( j, A 2 cr, A 3 a)dt (4) 

where 1 is the 2x2 unit matrix and o ee io 2 . The unimodularity condition on T enforces 
X\X 2 X 3 = 1. The black branes are parametrised by three charge parameters qj and take 
the form 



ds i = _#-2/3 fdt 2 + ^dr 2 + r 2 H l '*dxl 
Xi = , Aj 



(5) 



We have defined 

tf,(r) = l + ^, H(r) = H 1 H 2 H 3 , f(r) = - + H ^ +) +r 2 H{r) (6) 

and / = 1,2,3 labels each 2x2 block. The horizon is at r = r + ; we use the coordinate 
freedom in r to fix r + = 1 and thus Ki = qj in the conventions of [18]. 

The thermodynamics of these backgrounds was originally studied in [19, 20] in the context 
of holography, with the study of spinning D3-branes, their counterparts in the full 10D theory, 
initiated in [21, 22] and extended in [23, 24]. 

2.1 Multiple black brane branches 

We work in the grand canonical ensemble throughout. The chemical potentials for this 
system can be expressed in terms of the charge parameters 

Vgj(l + gl)(l + g 2 )(l + g3) 

= l + q~i ' (7) 
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and are conjugate to the three R-charges. In this ensemble the relevant thermodynamic 

n 

voi 3 ' 



potential is the Gibbs potential, Q, with density uj = We employ a dimensionless 



expression uj: 



(l + gi)(l + g 2 )(l + g 3 ) 16ttG 5 u; 

(*J — TT, 1 = 



16ttG 5 ' (mi + /4 + m§) 2 ' 

Similarly we employ a dimensionless version of the temperature, T: 

rp = 2 + gl + 92 + Q3 ~ gig2?3 rf = T /g\ 

27r v /(l + g 1 )(l + g 2 )(l + g 3 )' V^i + ^l + ^l' 

It is well known that this black brane family is not thermodynamically stable everywhere. 
For a stable thermodynamic equilibrium, the theory should lie at a local minimum of the 
thermodynamic potential, at fixed temperature and chemical potential. Constructing the 
Hessian matrix for variations of the energy density e, 

d 2 e 

n = asai, < 10 > 

with respect to the entropy density s and the R-charge densities pi following [18], a necessary 
condition for local thermodynamic stability is given by 

det H oc 2 - q x - q 2 - q 3 + gig 2 <?3 > 0. (11) 

However note that this is not a sufficient condition for local thermodynamic stability as 
there may be an even number of negative eigenvalues of H. 3 At a fixed set of chemical 
potentials, (/ii, /12, ^3), we find that there is a minimum temperature for the existence of a 
thermodynamically stable dominant branch of solutions. 

Without loss of generality we focus on fij > at which we only consider those solutions 
with qj > 0. This ensures that the solutions are regular. For this ensemble we fix two 
dimensionless ratios of the chemical potentials, Ci = ^ an d C2 = j^, or any function thereof, 
in addition to T. In general there are four black brane solutions at any grand canonical 
coordinate, (d,(2,T). To explore the phases of the system we adopt the coordinates 

/77 / 9 9 9 5 P*2, r— I q 9 9 \ / 

V6VM! + M2 + /i3 v2V^! + /^2 + ^3 



3 We wish to thank Mukund Rangamani for bringing this point to our attention. 
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which are dimensionless quantities and depend only on £ 12 . 

At a representative fixed (pi,P2) we show the free energy as a function of the temperature 
in the left panel of figure 1, illustrating that the dominant line of black branes are only 
stable above a threshold temperature. This is the threshold for thermodynamic instability 
as described by the condition (11), as well as the positivity of the eigenvalues of H. In 
appendix A we show that there is always one dominant branch of black branes with a 
minimum temperature given by the threshold for thermodynamic instability 



f f 




Fig. 1: Left: Gibbs free energy density for the various black brane branches that exist at fixed 
(PiiP2) = (—0.2,-0.1), defined in (12). Solutions that are thermodynamically unstable 
are shown in blue. A locally thermodynamically stable branch has a minimum temperature 
marked by the dashed line. Right: Fixed pi = P2 = 0. The equal-charge black brane is only 
thermodynamically stable above the temperature indicated by the dashed line, whilst below it 
is both thermodynamically unstable and subdominant in the grand canonical ensemble. 

It is worth noting that in some cases there is a branch of black branes which satisfy 
the condition (11), but which are not thermodynamically stable due to an even number of 
negative eigenvalues of "H. For example, the three equal chemical potential case, where pi = 
P2 = 0, is shown in the right panel of figure 1. The dominant black brane above the critical 
temperature T ~ 0.09 is locally thermodynamically stable and has equal charges q\ = g 2 = Q3 
and is nothing but the Reissner-Nordstrom-ArfS's (RN) black brane. Below this critical 
temperature % picks up two equal negative eigenvalues with eigenvectors 6(s, pi, p 2 , p 3 ) = 
(0, —1, 0, 1), (0, —1, 1, 0). These eigenvectors for this thermodynamic instability are naturally 
associated with a Gubser-Mitra dynamical instability [25, 26], as well as the generation of 
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neutral zero modes at the threshold of stability which we discuss further in section 5. 

On a different tack, the free energy analysis presented here reveals a worrisome issue 
associated with the RN solution, as considered in this top-down approach. Working in the 
further consistent truncation of our theory corresponding to three equal charge parameters, 
RN becomes the only black brane solution and therefore appears to be thermodynamically 
dominant. However, below the critical temperature for the thermodynamic instability dis- 
cussed above, T ~ 0.09, RN is not thermodynamically dominant in our more complete 
theory. 

To understand the coordinates (12) better, we project the condition fij > into (p\,P2) 
space. This region is formed by the intersection of three ellipses, with equations 

3p 2 x +£ = 1, 3p? + hp\ ± 2VSpm = 2 (13) 

and is shown in figure 2. Within this region we find there is always a stable branch that 
dominates the ensemble above some critical temperature, as in the explicit examples shown 
in figure 1. The minimum temperature for this dominant stable branch as a function of p\ 
and P2 is illustrated in the right panel of figure 2. 

3 Fluctuations about the normal phase 

We now detail the fluctuations of the 20 scalars and 15 gauge fields about the background 
normal phase solution given in section 2. We parameterise deviations from the scalar field 
matrix To (3) by the fluctuation matrix x( r )> an d deviations about the gauge field matrix 
Aq (4) with radial electric fluctuations, a(r): 

T = T e x , trx = ; A = (A + a)dt (14) 

where the trace condition enforces unimodularity of T at all orders in the fluctuation. The 
action for the linearised fluctuations can be found in appendix B. Since the background is 
block diagonal, the fluctuations in different 2x2 blocks do not couple. 

In general the fluctuations \ and a source metric fluctuations. For general background 
{qi} these may be switched off by enforcing the conditions 

tr xii = tr van = (15) 
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Pi Pi 

Fig. 2: Left: Region where black brane backgrounds with [ij > exist. Right: Critical temperatures 
for the thermodynamic instability — i.e. the minimum temperature, T, for which the dominant 
branch is thermodynamically stable. 

where the trace is performed within the 2x2 block. In detail these conditions are so restrictive 
that we have an = and a\n + Xn 17 — 0, leaving the only non-zero diagonal block 
fluctuation as o\n — Xn a i which as we shall see shortly is a charged mode. For special cases 
of the background the conditions relax: for two or three equal charges we find 

^ tr xn = tr aan = ( 16 ) 
i i 

where the sum is performed over equal qi blocks, with any remaining diagonal blocks con- 
strained by (15). 
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3.1 Diagonal fluctuations 

Without loss of generality we focus on fluctuations living in the (1, 1) block, 

aan \ 
(17) 
/ 

where xu is a symmetric 2x2 matrix and du is proportional to the 2x2 identity matrix. 
To satisfy the backreaction and unimodularity constraints away from equal charges we must 
impose tr^ii = 0. 

The equations of motion for these fluctuations may be partially diagonalised by moving 
to the new variables 

Neutral: <f>^ = axn + Xu ") a ii — &&11 + &ii&, (18) 
Charged: = a%\\ — Xno", oq^ — a&n — ana — 0. (19) 

The gauge potential for the (— ) mode vanishes due to the antisymmetry of aan. Note that 
in what follows we have specialised to radial dependence in deriving these equations, despite 
their being presented in partially covariantised form. 

3.1.1 Charged equations of motion 

The equation of motion for the charged scalar <f>^ is 

□4? - 4 (Xi (Xi - X 2 - X 3 ) + Al) <j>$ = 0. (20) 

Various special cases of this mode were considered nonlinearly in [5] for the one, two and 
three equal-charge background; for these cases we have numerically verified the critical tem- 
peratures quoted. The unimodularity and backreaction constraint tvxn = can be imposed 
without trivialising this mode. 
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3.1.2 Neutral equations of motion 

The equations of motion for the neutral scalar <j)^j together with the gauge field fluctuation 
are 




xf 1 x\ 



- AX^X, -X 2 - X 3 )0St } - + = (21) 

(22) 



where we have defined h = daffi . 

If the mode is to be considered alone as in (17), we must impose tvxu = for which we 
have = by (18). In order to satisfy the backreaction and unimodularity constraints 
without trivialising these modes, we must therefore include other diagonal blocks, but this 
can only achieved if two or more of the background charges are equal. Such a sum of modes 
at two equal charges is closely related to the Gubser-Mitra instability [25, 26] — a connection 
that we explore further in section 5. 

3.2 Off-diagonal fluctuations 

Without loss of generality we focus on fluctuations living in the (1, 2) block (and (2, 1) by 
symmetry) , 

/ X12 \ / aa 12 \ 

X = \ X T i2 , a = [ aj 2 a (23) 
\000/ \0 00/ 

for which the condition tr% = and backreaction constraints are automatically satisfied. 

As for the diagonal fluctuations, the equations of motion for these fluctuations may be 

partially diagonalised by moving to the new variables 

(+) mode: <f$ = a\\2 + XW, = a °<-n + (24) 

(-) mode: 0^ = a Xi2 - Xi2&, «i2 = ™i2 ~ (25) 

where in contrast to the diagonal modes, here both the (+) mode and (— ) mode are charged 
and couple to gauge field fluctuations in general. Dropping the indices and defining 

A ± = l(A lT A 2 ), F± = ^±^, h ± =da ± (26) 
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we obtain the coupled equations of motion 



□0± - ( - 2{X X + X 2 )X 3 + -F£ 

1 2 \ 



-g 



X1X2 



(2ht-{X l +X 2 )F^(j )± ) 



{X 1 -X 2 f tt , (X 1 + X 2 )(X 1 -X 2 ) 



X\X 2 



-g""a ± + 



When qi = q 2 we have that (f)\ 2 satisfies the same equation of motion as the charged 



X\X 2 



0. 



(27) 



(2f 



mode 4>\_i (20), and <f>^ } satisfies the same equation of motion as the neutral mode <f>^_ 
(21). 



a+) 



(+) 



3.3 Marginal modes 

Collecting all the modes above we have 

Type 



Charged diagonal 
Off diagonal (— ) 
Neutral diagonal 
Off diagonal (+) 



Scalar field See section 



*,(-) A-) A-) 
'u > ^22 > ^33 

6 {+) d {+) 
J n > V22 7 ^33 

^7^7^ 



4.1 
4.2 

5 
5 



Note that the equation of motion for the I 7^ J mode (j)fj reduces to the equation of motion 
for the (f>f^ when the background charges on the legs coincide: qj = qj. 

We seek time-independent solutions that are regular at the horizon and normalisable at 
the AdS boundary for each of the fluctuations detailed in the table above. Fluctuations of 
the gauge field at the boundary are required to vanish so that the chemical potentials are 
not perturbed. Such marginal modes indicate new branches of superconducting and normal 
black brane solutions. Developing series solutions we find the near-horizon behaviour (here 
the horizon position has been fixed to r + = 1) 



a 



% + . . . + log(r 
a 3 (r-l) + ... 



l)(a 2 + ...) 



(29) 
(30) 
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where the parameters a n are undetermined in the expansion and the ellipses denote deter- 
mined higher order terms. For regularity we set a 2 = 0, and since the system of equations 
is linear we are free to scale a± — 1. Additionally, for the diagonal modes we have fixed a 
shift symmetry in a such that a(r + ) = 0. The scalar fluctuations in the UV vacuum all have 
mass m 2 = —4, which saturates the Breitenlohner-Freedman bound there. The boundary 
series takes the form 



with data b n . Normalisability demands that we set b 2 = and we must also impose 63 = 
to preserve our choice of chemical potentials. 



In this section we consider fluctuations of the fields labelled by (— ). Before considering 
general ratios of the chemical potentials and general pi,P2 coordinates it is instructive to 
study a slice of the allowed region (figure 2) at /ii = fi 2 , for which p 2 = 0. We find marginal 
modes with critical temperatures illustrated in figure 3. 

On this slice we can consider the effect of an imbalance on the mode 4> 3 ~^\ which is 
charged under A 3 and does not couple directly to any other gauge field. At p 2 = 0, p\ 
becomes a function of the imbalance for this mode, £1,2 = The imbalance £1,2 = at 
Pi = —a/2/3, £1,2 = 1 at px — and £1,2 — > 00 as p\ — > a/1/3. As shown in figure 3 the 
critical temperature as a function of this imbalance meets the temperature of thermodynamic 
instability at sufficiently high £12. Note however that before this value of (1^ is reached we 
see superconducting branches for the modes 0i 1 \</4 2 ' ) and (f)^ , which are charged under 



Extending the analysis to the full allowed region of (^1,^2) space, we find the modes with 
the highest critical temperature as labelled in figure 4. 



(f) = r 2 (6i + ...) + r 2 logr(6 2 + • • •) 
a = b 3 + ...+r- 2 (&4 + •••) 



(31) 
(32) 



4 Sector (— ) 



.4 



12 



0.5 



0.4 



0.3 



0.2 



0.1 



0.0 



-0.8 



-0.6 



-0.4 



-0.2 



0.0 



0.2 



Pi 



0.4 



Fig. 3: A /ii = fi2 slice (/i3 ^ 0) as a warm-up to the full (pi,P2) space of normal phase solutions. 

The light blue region is thermodynamically unstable. The lines indicate critical temperatures 
of green: , red: 4>vs ■> rfhz < blue: 4>ii\ 4^22 ^ 4*12 • Finally we show the line of neutral 
scalar zero modes (see section 5) in black: 



J (-) 6 { - ] 6 { - ] 
'11 ' V22 ' Y12 

6 {+) 6 {+) 
'11 ' ^22 ' ^12 



4.1 4\4 ) :4 ) 

The critical temperatures for 33 ^ are plotted in figure 5. It is important to emphasise that 
the region where (f)^ is unstable is not restricted to the segment indicated in figure 4 - 
there is a slight overlap of <f)[i \ 4^ an d 4>2a near lines of equal chemical potential. 

4.2 <f$,<t$,t£ 

We find that marginal modes exist in a window of allowed p\ and P2 values. This region and 
their critical temperatures are plotted in figure 6. Note that the critical temperatures all 
lie below those found in section 4.1, except along the locus of equal chemical potentials on 
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Fig. 4: 




lodes which give the highest critical temperature. The diagonal modes 0^ , 4>22 > 033^ 
give the highest critical temperature in the blue regions as shown. In addition, along the red 
lines the modes 0^ , 0i 3 \ 023^ linearly coincide with the diagonal modes and have the same 
critical temperature. Dots indicate special cases of the chemical potentials, (^1,^2,^3), with 
all equal at the origin. 



the two legs; for instance along p 2 = the critical temperatures for (f> l2 and an d 022 
coincide. 

5 Sector (+) and Gubser-Mitra instabilities 

In this section we consider fluctuations of the fields labelled by (+). The fields 0!q , 2 2 , 033 
are neutral, while 4>^\ 0i3~\ 023 become neutral when two chemical potentials coincide. We 
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P2 0.0 




Fig. 5: Critical temperatures for the mode 33 \ The modes 4h.li ^22 can ^ e obtained by permuting 
the labels and appear predominantly in the regions as labelled in figure 4. 



find marginal modes along the loci 





6 i+) 




/'I 


= M2 > M3 


it', 


6 i+) 
V33 > 




Ml 


= M3 > M2 


(+) 

22 > 


< T33 1 


& 


M2 


= M3 > Ml 



(33) 



where all fluctuations are neutral. The backreaction and unimodularity constraints for the 
diagonal modes can be satisfied by turning on two modes, e.g. <f)^ = —<Pxi- 

It is perhaps surprising that unlike the (— ) modes of section 4, these marginal modes 
exist only for one dimensional regions of the {p\,P2} plane. However, they do not indicate a 
true instability of the black brane background. Instead, as we shall demonstrate, they have a 
natural interpretation as zero modes which move in the family of background black branes. 
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Fig. 6: Region where marginal modes for <p\ 2 exist, with critical temperatures indicated by the 
contours. Regions for the other modes </4 3 ^ and <^ 3 can be found by permuting the indices. 
The blue shaded region is the region shown in figure 2. 



The existence of these zero modes is closely related to the threshold of thermodynamic 
instability; indeed, the critical temperatures are given precisely by the lowest temperature 
at which the normal phase is locally thermodynamically stable. 

We can analytically construct this class of marginal modes in the grand canonical en- 
semble by perturbing the charges and the horizon position whilst ensuring that the chemical 
potential and temperature remain fixed. A 4x4 matrix for this variation is 

M = 'p&r (34) 

and, if the fluctuations are to be marginal modes in the grand canonical ensemble, we 
require det Ai = 0, which is achieved for background values of qi at the threshold of the 
thermodynamic instability given by det% = 0. We emphasise that this applies to modes 
which include gravitational fluctuations at any values of the chemical potentials and is easily 
derived in this way. For example, when q% = > we have a surface of zero modes 
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generated by perturbations to (5), Hj — > Hj + 6Hj, f — > f + 5f, where 



SHi = r^^rv^, i = l,2 (35) 

(qi - 1) r 2 

r tt Jl + 92-2 1 

SE > = -fe-Dfe-i)? (36) 

= 2+ + a -2) 

Specialising to the non-gravitational modes of interest here which exist only along loci of the 
chemical potentials, we obtain 4 



8X\ 5X2 



X 1 X 2 1 + r 



2 ■ 



SX 3 = (3* 



SA 1 = -SA 2 = ^±M^i) , Ms = (39) 

(1 + r 2 ) 

at temperature T = — , 1 =. This zero mode corresponds to a redistribution of the 

7r v /2(l+4g 3 +'?I) 

black brane charges when two of them are equal. These modes satisfy equations (21) and 
(22) and consequently it is precisely this mode which gives rise to the black line in figure 3. 
It is also the zero mode generated by the thermodynamically unstable directions discussed in 
section 2.1. We expect that there is also a true instability associated with these fluctuations 
[25, 26] but it is not captured by a zero mode analysis. We also anticipate a true instability 
away from the lines of equal chemical potentials, provided gravitational fluctuations are 
included. 



6 Discussion 

We have found new branches of superconducting black branes in Af = 8 £0(6) gauged 
supergravity, coming from a truncation involving 20 scalar fields and 15 gauge fields. Su- 
perconductivity persists to any values of the chemical potentials conjugate to the three 

4 The no-backreaction condition imposes that any two of the background charges are equal; let us here 
choose qi = q%. Subsequent conditions on the fluctuations themselves become 5q\ + Sq2 = and Sq^ = 0. 
Further imposing that the chemical potential is not perturbed we have q2 = f . 
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R-charge densities of the background black brane. The critical temperatures at which these 
black brane solutions arise is always higher than the temperature at which the three-charge 
background becomes thermodynamically unstable. 

It is known that in some special cases — such as two and three non-zero equal chemical 
potentials — one has 'retrograde condensation', i.e. the branch of superconducting black 
branes is thermodynamically subdominant and the condensate exists only above a critical 
temperature [5] (see also [27, 3] in 4D). 5 The other previously studied case of this system is 
that of one non-zero chemical potential, where the superconducting black brane is thermo- 
dynamically dominant [5]. Here, in all known cases of retrograde condensation, we find that 
there is an additional marginal mode with the same critical temperature (see figure 4). We 
leave a nonlinear investigation of this branch, and whether it is dominant, to future work. 

For a given mode we may interpret dialling the chemical potentials as introducing an 
imbalanced mixture as initiated in [12, 13]. The critical temperature surface cannot be 
extended to arbitrarily low temperatures as the imbalance is dialled, due the presence of 
a thermodynamic instability in the normal phase at finite temperature. At a given set 
of chemical potentials, different modes in the consistent truncation experience a different 
effective imbalance. Interestingly, we find that before a superconducting instability vanishes 
as a function of its imbalance, another appears to take its place. 

The difference of the chemical potentials sets a scale for the system. It would be inter- 
esting to see if this scale sets a value of spatial momentum in this truncation resulting in 
an inhomogeneous, FFLO-like phase [29]. 6 It would also be interesting to construct all of 
the superconducting branches found nonlinearly, away from the critical temperature. This 
would allow the direct comparison of the free energy with the normal phase, and consequently 
the identification of the true phase boundary for the spontaneous breaking of a [/(l)-bulk 
symmetry. 

For many of the linear modes studied in this paper, labelled by (+), we have found 
only neutral marginal modes, whose role is to move one infinitesimally within the family 
of background black branes. These modes occur precisely at the threshold temperature for 

5 For systems exhibiting retrograde condensation, the expectation is that superconducting ground states 
do not exist in the planar theory, but can exist when the boundary has spherical topology [28]. 
6 For examples of spatially modulated phases in a truncation of IIB supergravity on S 5 see [30] . 
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thermodynamic instability. As such the marginal modes do not generate an actual instability 
of the black brane. We have shown that it is possible to populate the entire minimum 
temperature surface (as a function of chemical potentials) with zero modes of the background 
black brane, which can be found analytically and involve gravitational fluctuations in general. 
One expects to find a true instability associated with the thermodynamic instability in these 
cases. 

Our work was motivated in part by the search for universal features in top-down holo- 
graphic superconductivity. We wish to highlight a basic problem we have encountered when 
using consistent truncations: a subdominant solution in a given theory can appear as the 
dominant solution in a consistent truncation of that theory. This is exemplified by the 
normal phase solutions employed in this paper. 

In the case of three equal chemical potentials the normal phase solution is simply RN at 
high temperatures. Whilst the RN branch exists at low temperatures too, it is no longer the 
dominant contribution to the ensemble. This has important consequences for the equal-g/ 
consistent truncation where RN is the only normal phase solution remaining. Evidently in 
the equal-g/ truncation RN appears as the dominant black brane branch at all temperatures, 
when in fact, at low temperatures in the full theory it is subdominant to a thermodynamically 
unstable black brane. 
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A Dominant branch 

Using q 3 as a parameter along the branch of solutions and writing sj = (1 + g 3 ) 2 — Aq 3 Q we 
find that qi , q 2 are given by 



Qi = 



1 + <?3 ~ Si 
1 + <?3 + «1 



Q2 = 



1 + <?3 ~ S 2 

l + q 3 + s 2 



(40) 
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There are therefore up to four solutions coming from the sign choices for si and s 2 . In these 
coordinates the thermodynamic potential density becomes 

(l + q 3 )(l + q 3 + Sl )(l + q 3 + s 2 ) 

4?I(i + Ci +C!) 2 1 } 

and so we see that at some desired grand canonical coordinate, (d,(2,T), the dominant 
configurations live on the s±,s 2 > branch, provided it exists there. The temperature is 
stationary at the threshold for thermodynamic stability: 



dT 

dq 3 



l-q 2 3 )(s 1 + s 2 ) + 2s 1 s 2 ^ 



CiCa 87rg 3 si W?s(l + Ci + CD 



det^ oc 2 ^-^ + s,) + 2s lS2 
(l + q 3 + Sl )(l + q 3 + s 2 ) 

and furthermore, on the S\,s 2 > branch the temperature is a minimum there: |^ > 0. 
This shows that there is a branch of black branes which only exists above the threshold 
temperature for thermodynamic stability, where it dominates the ensemble. 

B Action for fluctuations 

In this appendix we present the action to quadratic order for fluctuations of (1) of the form 
(14) that do not source metric fluctuations: 



l — I d 5 xy/^g~ tr (L x + L a + L int ) (44) 
i"G 5 J 



with 



£2 — 1 , 

i07T 



L x = - 9 ^x' 2 - 9 ^((A oX r-Al X 2 ) + 9 ^((T^A^ (45) 



-9 2 c {t 2 ,-\^T )T^ x 2 -9 2 (ToX) 2 



L a = ~ 9 ^( a 2 -T^aT a) + 9 -^f(T^ar (46) 
L^t = - 9 ^T 1 aT (A oX -xAo)- 9 ^{xT 1 A' + T 1 A' x)T - 1 a' (47) 
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